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The determination of dark matter constraints from liquid xenon direct detection experiments de- 
pends upon the amount of scintillation light produced by nuclear recoils in the detector, a quan- 
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Figure 1: Various measurements of the scintillation efficiency factor «S? e ff ■ Figure courtesy of M. Schumann. 

1. Introduction 

The DAMA annual modulation [1,2] and the excess low-energy events in CoGeNT [3] can be 
explained by recoils from a light Weakly Interacting Massive Particle (WIMP); see e.g. Ref. [4]. 
The dark matter interpretation of these signals can potentially be confirmed or refuted by exper- 
iments such as XENON10 [5] and XENON100 [6, 7], but the ability of these experiments to 
probe light WIMPs is dependent on their sensitivity to low energy recoil events in their detec- 
tors. The calibration of these detectors' energy scales (and thus knowledge of their sensitivity to 
low-energy events) is dependent upon the scintillation efficiency factor -£? e ff. We examine here 
how the XENON constraints depend upon various models for the poorly known ^f e ff - 

Some of the work presented here is described in Ref. [8], although the analysis of the XENON10 
results has been extended down to lower recoil energies with a more careful treatment [9]. 

2. The =£? e ff Scintillation Efficiency Factor 

Recoils in the XENON detectors produce two scintillation signals, referred to as 51 and 52, 
resulting respectively from prompt photons and ionization produced by the collision of a WIMP 
with a xenon nucleus. Interpretation of the XENON results requires the ability to reliably recon- 
struct the nuclear recoil energy E nr from the observed 51 signal. Calibration of the nuclear recoil 
energy dependence of 51 often involves gauging the detector's response to electron recoils at higher 
energies; parts of the detector's response (e.g. the fraction of scintillation photons that yield pho- 
toelectrons (PE) in the photodetectors) are more easily determined in this case than with nuclear 
recoils at lower energies. Taking SI to be normalized to the number of PE, SI and E nr are related 
by an equation involving the higher energy electron recoil calibrations: 

(51) = (S nr /S ee )£> eS (E nr )LyE nr . (2.1) 

Here, L y is the light yield in PE/keVee for 122 keV /-rays. ^£^{E nr ) is the scintillation efficiency 
of nuclear recoils relative to 122 keV y-rays in zero electric field; this factor is a function of the 
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Figure 2: (left) Assumed models for Jz? e ff, based on the Manzur et al. measurements [11] with three different 
extrapolations below E nr = 3.9 keV: constant (blue), linearly falling (red), and zero (green). Solid curves 
indicate fiducial cases, while lighter regions indicate la uncertainties, (right) The expected average SI 
signals in XENON10 for these Jz? e ff models. XENON100 exhibits a similar (but not identical) behavior. 



nuclear recoil energy. Since there is an applied electric field in the experiment, which reduces the 
scintillation yield by quickly removing charged particles from the original interaction region, two 
additional factors must be taken into account: S ee and S nr are the suppression in the scintillation 
yield for electronic and nuclear recoils, respectively, due to the presence of the electric field in the 
detector volume. The quantities S ee , S nr , and L y are detector dependent; Jz? e ff is not. Equation (2.1) 
describes the average SI signal ((SI)); the actual observed SI signals exhibits some random fluc- 
tuations about this average as described in Section 3. 

A variety of -Sf e ff measurements and estimates are shown in Figure 1 . These measurements 
are plagued by systematic issues and give conflicting estimates of -Sf e ff at low recoil energies; see 
Ref. [10] for a discussion regarding the various issues in these Jz? e ff measurements. There are two 
issues here: (1) Which of the Jz? e ff measurements should be used as a basis for analyzing direct 
detection results? and (2) Measurements of Jz? e ff have only been made at energies above some 
minimum; what is the behavior of J2? e ff at low energies, where no measurements have as yet been 
made? We do not address the first issue (which set of measurements to use) and simply choose 
the set of fixed-energy measurements that will give the most conservative constraints: the Manzur 
et al. measurements [11, 12] (black points). For the second issue (the behavior of J§f e ff at energies 
below the existing measurements), we consider three extrapolations of Jz? e ff below recoil energies 
of 3.9 keV (the lowest Manzur measurement): (1) Jz? e ff is constant, (2) Jz? e ff is linearly falling with 
recoil energy, and (3) ^? e ff is strictly zero. The XENON estimates of Jz? e ff have been suggestive 
of the first (constant) case, while the Manzur measurements are more suggestive of the second 
(falling) case. The third case leads to the most conservative possible constraints by simply ignoring 
any contributions from recoils with energies below 3.9 keV and is not necessarily a realistic model 
of the low-energy -S? e ff behavior. In addition to the dependence of the XENON constraints on these 
low-energy extrapolations, we consider how the constraints vary within the la uncertainties in the 
Manzur measurements. Our ^f e ff models are shown in the left panel of Figure 2, while the right 
panel of the figure shows the corresponding average 51 signals in the XENON10 detector. 
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Figure 3: (left) Monte carloed events in the XENON 10 detector for nuclear recoils with energies that cor- 
respond to average SI signals {(SI)) of 1.0 and 10.0 PE (black dots). Events that fail to be identified via the 
SI peak finding process are shown in gray. Events that are identified but fall either below the S2 threshold 
(black curve) or outside the nuclear recoil band cut (red curves) are shown in green. Events that are identified 
and pass all cuts (signal events) are shown in purple, (right) The fraction of events at each recoil energy that 
will be identified and pass all cuts in XENON 10 (2 < 51 < 75 PE). The three colored curves at low energies 
correspond to the three low-energy «5? e ff extrapolations shown in Figure 2. 



3. Detection in XENON10/100 



The standard analysis of XENON 10 and XENON 100 results uses the two scintillation signals 
SI and S2. Due to various physical processes and detector limitations, the observed signals exhibit 
random fluctuations about the average expected values at a given nuclear recoil energy. The most 
significant (but not only significant) contribution to the fluctuations is the Poisson fluctuations in 
the small number of PE's produced by the prompt scintillation (only ~10% of the prompt photons 
yield a PE), which leads to a poor energy resolution at low recoil energies. We show in the left 
panel of Figure 3 a monte carlo example of the scattering of the observed SI and S2 scintillation 
signals in XENON 10 for recoil energies yielding an average SI of 1.0 and 10.0 PE. We also show 
the two main cuts — an S2 threshold and the nuclear recoil band cut — and take into account the two- 
fold PMT detection efficiency (also known as the SI peak finding efficiency factor r/^i ). We use 
the low-threshold (SI > 2.0 PE) XENON10 results of Ref. [5] in our analysis (see also Ref. [13] 
for a detailed description of the XENON10 detector). We stress that this is a threshold in SI, not 
in (SI): as shown in Figure 3, a significant number of events produced at a recoil energy yielding 
(SI) = 1.0 PE will have an observed SI > 2.0 PE and pass all cuts, so low-energy recoils can still 
contribute to the signal within the XENON10 analysis region. To determine the contribution of 
these low-energy recoils in XENON10, we monte carlo a large number of events at each recoil 
energy and determine the fraction of events that will be both observed and pass all cuts [9] ; the 
results are shown in the right panel of Figure 3. See Ref. [14] for further discussion regarding the 
fluctuations of the signals in XENON10 and issues in determining the detection efficiencies. 

XENON100 has a lower background than XENON10, but has a higher threshold of SI > 4 PE 
[6, 7]. As the full details of the various data cuts in XENON100 are not yet available, we use for 
this experiment a simpler model of the efficiency and impose a (SI) > 1.0 PE cutoff as described 
in Ref. [8]. 
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Figure 4: XENON 10 (green) and XENON 100 (purple) 90% C.L. exclusion constraints for a constant Jz? eff 
at recoil energies below 3.9 keVnr. The solid curves are the constraints using the central values of -£? e ff; 
dashed curves and lighter filled regions indicate how these 90% constraints vary with the 1 a uncertainties 
in ^ e ff- Overlapping XENON10 and XENON100 regions are shown in blue. Also shown are the DAMA 
modulation 90%/3o75(7-compatible regions (gray contours/region) and the CoGeNT 90%-compatible region 
(pink contour/region over 7-12 GeV). 
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Figure 5: Same as Figure 4, but taking Jz? e ff to fall linearly to zero for recoil energies below 3.9 keVnr. 

4. Results and Discussion 

We now examine how the XENON 10/1 00 WIMP mass & cross-section constraints are affected 
by the ^f e ff model. We consider only spin-independent elastic scattering and assume an isothermal 
halo model as described in Ref. [8]. The results discussed below apply only to this case; spin- 
dependent couplings and alternate halo models may affect the compatibility of various experimental 
results. 

Our main results are shown in Figs. 4-6, corresponding to the three cases for the behavior 
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Figure 6: Same as Figure 4, but taking Jz? e ff to be zero for recoil energies below 3.9 keVnr. 

of -£? e ff at low recoil energies. The 90% confidence level (CL) exclusion limits for the fiducial 
(central) Jz? e ff models are shown as solid curves for XENONIO (green) and XENON100 (purple); 
lighter filled regions correspond to how these 90% CL exclusion curves vary with the 1 a uncer- 
tainty bands in the J2? e ff models. For comparison, we also show the WIMP parameters compatible 
with the DAMA modulation [1,8] (gray contours/regions, corresponding to compatibility within 
the 5(7, 3(7, and 90% CLs) and the region suggested by CoGeNT [3] (pink contour/region, corre- 
sponding to compatibility within the 90% CL). 

The XENON100 constraints are nearly identical for all three J2? e ff models. For the fiducial 
cases, XENON100 excludes all of the DAMA 3a region, indicating signficant incompatibility 
between these two experimental results for this dark matter model. XENON 100 excludes only the 
9-12 GeV WIMP mass part of the CoGeNT region, allowing for 7-9 GeV WIMPs. When the la 
variations in the Manzur data are considered, however, XENON 100 can exclude as little as only the 
DAMA 90% CL region and almost none of the CoGeNT 90% region to as much as nearly all of the 
DAMA 5a CL and CoGeNT 90% CL regions. The similarity of the XENON100 constraints for 
the different Jz? e ff models is due to the 4 PE threshold: at recoil energies below 3.9 keV, where the 
models differ, ~ 1 or less PE are expected on average, so recoils at these energies are unimportant. 

Due to the lower threshold of 2 PE, the XENONIO constraints are more dependent on the 
Jz? e g" model. The XENONIO constraints for the three fiducial -S? e ff cases vary significantly, with 
the low WIMP mass cut-off varying between ~4 GeV and ~7 GeV. In all three cases, though, the 
XENONIO constraints exclude the entire CoGeNT 90% CL and DAMA 3a CL regions. Though 
the average 51 signal is below threshold for recoil energies below 3.9 keV (where the Jz? e ff models 
differ), the random fluctuations in the SI and S2 signals lead to some events at lower recoil energies 
being observed above the threshold. For the light WIMPs in the region of interest for DAMA and 
CoGeNT, the number of low-energy recoils is very large and, even if the fluctuations are small, a 
significant number of events should still be observed above threshold. If the la bands on Jz? e ff are 
taken into account, XENONIO still excludes all of the DAMA 3a region and all but a narrow band 
of the CoGeNT region over WIMP masses of 7-9 GeV. 
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While the ability of the XENON dual-signal (SI and 52) analyses to constrain light WIMPs 
is particularly sensitive to systematic issues such as the behavior of -S? e ff , we note that this type 
of analysis is optimized more for heavier WIMPs. An analysis based on only the 52 signal, such 
as presented by P. Sorensen in this conference [15], provides much greater sensitivity to light 
WIMPs. However, we have shown here that the standard dual-signal analysis can still provide 
stringent constraints on low mass WIMPs even for conservative models of Jz? e ff - 
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